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INTRODUCTION 


Recent  analytical  investigations  suggest  that  a  reduction  in  wall  bounded  turbulent  drag  can  be 
achieved  by  the  action  of  a  transverse  traveling  wave  via  surface  deflections.  The  reduction  in 
drag  can  be  up  to  40%  which  could  eventually  lead  to  significant  improvements  in  aircraft 
efficiency.  This  work  has  mostly  been  verified  analytically  since  the  actuator  parameters  (forces, 
displacement,  and  frequency  range)  required  for  this  traveling  wave  are  beyond  standard 
actuator  capabilities. 


This  effort  pursued  the  investigation  and  development  of  structural  waveguides  that  have 
potential  to  meet  these  specifications  and  experimentally  explored  the  ability  for  drag  reduction. 
A  structural  waveguide  can  be  theoretically  produced  in  any  structure  by  controlling  the 
boundary  conditions  to  make  the  structure  look  “continuous”.  This  can  be  accomplished  by  a 
combination  of  intelligent  structural  design  coupled  with  piezoelectric,  piezopolymer,  electro- 
expulsive,  or  other  embedded  smart  material  actuators.  Circular  structures  are  well-suited  for 
propagating  waves  since  they  are  inherently  continuous.  This  leads  directly  to  application  on 
aircraft  fuselages,  where  the  propagating  wave  in  the  circular  structure  will  be  inherently  parallel 
to  the  flow.  Smart  material  actuators  using  phased  array  drive  concepts  can  be  applied  to  the 
structure  to  excite  the  traveling  waves.  The  structural  waveguide  concepts  were  analytically  and 
experimentally  investigated.  This  was  a  one  year  proof-of  concept  preliminary  study. 

The  design  and  fabrication  of  the  structural  waveguides  was  completed  successfully  and 
structural  characterization  of  performance  confirmed  the  expected  results.  A  summary  of  these 
results  will  be  provided  in  the  flowing  sections.  At  the  end  of  the  project,  preliminary 
aerodynamic  validation  of  the  concepts  was  pursued  on  a  cylindrical  structure  instrumented  with 
a  full  array  of  flow  visualization,  vibration  and  drag  measurement  equipment.  Upon  completion 
of  the  performance  period  of  the  project  and  during  the  authorship  of  the  present  report,  the 
aerodynamic  experiments  were  still  in  process.  Further  experiments  analysis  and  detailed 
dissemination  of  the  aerodynamic  testing  results  will  be  presented  in  future/upcomming 
conference  proceedings  and  archival  journal  publications. 

This  was  a  multidisciplinary  effort  between  faculty  of  the  Mechanical  engineering  Department  at 
Virginia  Tech.  Drs  James  Carneal  and  Marty  Johnson  of  the  Vibrations  and  Acoustics  Lab 
provided  expertise  in  structural  dynamics,  phased  array  actuators/sensors,  and  novel  actuators 
for  active  noise  and  vibration  control,  and  Dr  Pavlos  Vlachos  provided  the  expertise  in  fluid 
dynamics  and  experimental  flow  diagnostics. 

The  research  presented  herein  is  original  and  risky.  The  preliminary  results  indicate  that  the 
proposed  concept  is  feasible  and  offer  scientific  and  practical  engineering  merit.  Continuation  of 
this  effort  has  the  potential  for  significant  scientific,  social  and  economic  impact  and  commercial 
applications.  Almost  all  wall  bounded  flows  and  fluid  machinery  are  subject  to  wall  turbulence. 
Efficient  control  of  the  character  of  the  flow  will  save  energy,  improve  performance,  and  open 
new  prospects  for  engineering  applications. _ 


BACKGROUND: 


In  the  following  section,  background  work  is  summarized  and  presented  in  two  main  areas: 
turbulence  control  theory  and  structural  waveguides.  However,  extensive  review  of  this  topics  is 
beyond  the  scope  of  this  document  as  well  as  review  of  topics  related  to  smart  material 
actuators,  and  phased  array  actuators. 

TURBULENCE  CONTROL  THEORY 

Numerous  flow  control  methodologies  have  been  employed  to  favorably  alter  flow  behavior. 
Applications  include  flow  control  of  separation,  separated  flows,  circulation,  as  well  as  lift 
increase  and  drag  reduction.  All  these  technologies,  upon  maturation,  can  potentially  have 
significant  impact  to  a  variety  of  vehicle  platforms  and  fluid  machinery.  Among  these,  turbulent 
boundary  layer  control  towards  viscous  drag  reduction  has  been  a  formidable  area  of  research. 
A  thorough  review  of  this  topic  is  beyond  the  scope  and  the  page  limitations  of  this  text.  We  will 
however  review  and  discuss  some  of  the  most  recent  and  significant  contributions  as  they 
pertain  to  the  effort  proposed  herein. 


The  development  of  modern  turbulence  control  mechanisms  and  the  corresponding  shear 
stress  reduction  depends  upon  the  observation  that  organized  structures  in  the  near-wall  region 
and  the  friction  drag  are  inevitably  linked  .  Two  key  elements  of  this  view  have  been  the 
discovery  of  low 


Figure  0:  Instantaneous  flow  of  wall 

speed  streaks  in  turbulent  boundary  ^ 


streaks.  Top:  no-control;  Bottom: 

layers  as  well  as  the  verification  of 


Theodorsen's  hypothesis  on  the 


Traveling  wave  excitation  corresponding 


to  1=1,  A  *=840  and  T*=50,  with  L  =840. 

existence  of  stable  three-dimensional  z  vortical  structures,  called 

.  .  2  (Du  &  Karniadakis,  1999) 

hairpins  .  More  specifically,  vortex  pairs  (loops)  and  horseshoe  or  hairpin-like 

vortices  have  been  recognized  as  the  dominant  flow  structures  in  the  near-wall  region  of 
turbulent  boundary  layers,  and  many  researchers  have  colorfully  described  boundary  layers  as 
a  “forest”  or  a  “jungle”  of  hairpins  . 


In  this  effort  we  are  interested  in  the 
mechanisms  of  modifying  the  above- 
described  structure  of  the  turbulent  boundary 
layer.  Karniadakis  &  Choi  provide  the  most 
up-to-date  and  extensive  review  on  flow 
control  of  turbulent  boundary  layers  via 
transverse  actuations.  Methods  for  near-wall 
turbulence  control  have  been  developed 
under  the  assumption  that  the  turbulence 
production  cycle  could  be  favorably  altered 
by  the  manipulation,  stabilization  or 
suppression  of  low-speed  streaks,  quasi- 
streamwise  vortices  or  the  hairpin-like 
structures  that  populate  the  near-wall  region. 
Small  grooves  or  riblets  mounted  on  the  wall 
surface  have  proved  to  be  effective  in 
partially  suppressing  turbulence  and  reducing 
the  drag  force  by  about  5%  to  10%  . 
Additional  promise  was  demonstrated  by  flow 
control  using  wall  oscillations  as  an 
alternative  effective  mechanism  for  drag 
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reduction  .  Drag  reduction  as  much  as  40% 

was  demonstrated  experimentally  by  Choi  et 

10 

al 


Despite  successful  demonstration  of  drag  reduction  methodologies,  the  fundamental 
understanding  of  the  governing  dynamics  remains  a  great  challenge.  Significant  contributions 
by  Jimenez  &  Pinelli  and  Schoppa  &  Hussain  have  shed  light  on  the  self-sustained  processes 
of  wall-turbulence  and  demonstrated  numerically  that  the  turbulence  regeneration  cycle  is 
exclusively  associated  with  the  near-wall  region  and  is  independent  of  the  outer  flow.  Schoppa 
&  Hussain  later  used  this  concept  to  propose  a  large-scale  drag  reduction  technique  employing 
vortex  generators  along  the  channel  centerline,  which  resulted  in  drag  reduction  on  the  order  of 
20-60%  .  Subsequently,  they  were  able  to  identify  a  new  mechanism  for  the  transient  growth  of 
streak  instabilities  .  These  contributions,  among  others,  provide  an  established  relationship 
between  streamwise  vorticity,  streak  instability  and  vortex-induced  viscous  drag.  Subsequent 
work  revealed  a  fundamental  mechanism  which  completely  eliminates  these  near-wall 
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streaks  .  Using  numerical  simulations  Du  &  Karniadakis  proposed  a  flow  control  method 
based  on  the  application  of  a  body  force  acting  within  the  viscous  sublayer  and  decaying 
exponentially  away  from  the  wall,  traversing  the  flow  in  the  spanwise  direction  in  the  form  of  a 
traveling  wave.  In  particular,  the  following  traveling  wave  force,  Fz,  is  introduced: 

y 

2^  2^  F  =  Ie  *sin( 


t)  Equation  1 

z 


^  T 

where  I  is  the  amplitude  of  excitation,  A  is  the  wavelength  (along  the  span),  and  T  is  the  time 

♦ 

period.  Here  A  =U*  A  /v,  where  U*  is  the  wall  shear  velocity  and  v  is  the  kinematic  viscosity  (“+" 
as  a  superscript  indicates  that  the  quantity  is  expressed  in  viscous  wall  units). 

For  channel  flow,  it  is  shown  in  Figure  0  (z  :  spanwise  direction)  that  the  uncontrolled  case 


exhibits  the  familiar  turbulence  structure  consisting  of  pairs  of  high-speed  and  low-speed 
streaks,  but  no  such  pairs  can  be  found  near  the  controlled  surface.  The  resulting  overall  drag 
reduction  from  this  simulation  was  in  excess  of  30%.  This  is  a  significant  finding,  in  view  of  the 
fact  that  streaks  and  streak  spacing  are  usually  difficult  to  alter,  even  in  cases  where  a  very 
large  amount  of  drag  reduction  has  been  obtained  as  in  drag  reduction  using  polymer  or  bubble 
injections.  All  previous  methods  were  unable  to  eliminate  the  near  wall  streaks,  whereas 
this  technique  accomplishes  exactly  that. 

Addressing  the  energetics  associated  with  this  turbulence  control  methodology,  it  has  been 
found  via  systematic  DNS  experimentation  that  the  specific  waveform,  i.e.  the  parameters  T,  A,  I,  are 
critical  in  obtaining  turbulence  suppression  and  corresponding  drag  reduction.  They 
hypothesized  that  it  is  the  total  amount  of  energy  input  that  is  critical  in  suppressing  turbulence, 
and  this  can  be  expressed  by  the  product  I  x  T .  This  quantity  represents  an  energy  invariant 
parameter  and  it  has  been  demonstrated  that  there  is  an  optimum  regime  of  energy  input, 
outside  which  the  control  has  either  no  effect  or  results  in  promoting  turbulence.  For  practical 
applications  this  is  very  important  as  it  suggests  that  the  flow  control  is  not  only  robust  but  is 
also  forgiving  in  the  sense  that  there  is  not  a  unique  input  condition  that  will  have  positive 
ramifications.  Any  product  of  I  and  T  that  satisfies  the  energy  requirement  should  be  effective. 
The  implications  of  this  are  that  this  aerodynamic  flow  control  method  is  open-loop, 
i.e.  the  control  system  does  not  need  estimates  of  the  plant,  nor  does  it  need  upstream 
sensors  to  correctly  phase  the  structural  waveguide  with  the  flow  perturbations.  It 
appears  that  this  optimum  regime  is  a  function  of  the  Reynolds  number  and  as  a  result  for 
optimum  control  the  penetration  length  needs  to  be  determined  from  the  reciprocal  of  the 
energy  invariant  product : 


A»  Const  x  (I  x  T+)-1 .  Equation  2 

From  Du  \  these  non-dimensional  parameters  can  be  applied  to  aerodynamic  turbulent  drag 
reduction,  resulting  in  wavelengths  on  the  order  of  20-200  mm,  and  pulsing  frequencies  on  the 
order  of  50-500  Flz. 

The  hypothesis  of  this  invariant  product  has  been  tested  extensively  numerically,  however  to 
date  no  experimental  proof  has  been  established.  There  is  very  little  information  available  for 
experiments  employing  an  electromagnetic  traveling  wave  force.  Karniadakis  and  Choi  state 
that  there  were  no  journal  papers  available  at  that  time;  however  they  report  on  a  few  ongoing 
preliminary  efforts.  Recently  Park  et  al  and  Breuer  et  al  reported  on  Lorentz  force  actuator 
development  and  drag  reduction  experiments  using  direct  force  balance  measurements.  This 
preliminary  study  was  limited  to  very  low  free-stream  velocities,  and  a  drag  reduction  on  the 
order  of  9%  was  found.  However,  the  mechanical  efficiency  was  approximately  10  which  is 
unacceptable  for  any  practical  applications.  Their  detailed  study  illustrated  several  problems 
related  to  the  implementation  of  a  Lorenz  force  actuation  scheme  such  as  corrosion,  bubble 
formation,  and  low  induced  velocities.  Reportedly,  the  efficiency  increases  for  high  Reynolds 
numbers  (Prof.  Karniadakis,  personal  communication);  however,  these  efforts  are  still  in 
progress  and  many  questions  remain  unanswered. 

Inspired  by  these  previous  efforts,  this  work  explores  an  alternative  methodology  for  turbulent 
boundary  layer  flow  control  that  overcomes  the  limitations  of  the  Lorenz  force  actuation.  A  flow 
disturbance  in  the  form  of  a  traveling  wave  force  generated  by  a  structural  waveguide  is 
used  to  manage  the  near  wall  turbulent  flow.  This  innovative  concept  overcomes  some  of 
the  shortcomings  of  the  electromagnetic  force  such  as  corrosion  or  bubble  formation,  and  more 


importantly  is  independent  of  the  conductivity  of  the  fluids,  thus  applicable  to  both  air  and  water 
flows.  Although  implementation  concepts  have  been  proposed  before  ,  no  published  results, 
experimental  or  numerical,  of  turbulent  boundary  layer  flow  control  using  these  concepts 
currently  exist.  Moreover,  the  mechanism  of  the  proposed  flow  control  is  based  on  modifying  or 
interrupting  the  self-sustaining  cycle  of  near-wall  turbulence.  Therefore,  in  contrast  to  previous 
successful  flow  control  approaches,  the  present  methodology  tackles  the  fundamental 
mechanisms  of  near  wall  turbulence. 

STRUCTURAL  WAVEGUIDE  CONCEPTS 

Analysis  of  a  structural  waveguide  is  a  fundamental  analysis  of  distributed  elastic  system 
flexural  vibration,  where  the  beam/plate  is  assumed  infinite  and  a  wave  is  free  to  propagate  in 
all  directions.  The  inclusion  of  finite  boundary  condition  to  this  model  will  result  in  modal 
behavior,  which  is  a  result  of  the  wave  bouncing  off  the  boundaries,  which  occurs  when  the 
propagating  wave  encounters  a  change  in  impedance.  Physically,  structural  waveguides  exist  in 
circular  structures  since  they  are  continuous  and  have  no  boundary  .  From  Fahy  ,  the  axial, 
tangential,  and  radial  displacements  (u,v,w)  of  the  wall  of  a  cylinder  vary  with  the  axial  position  z 
and  the  azimuth  angle  0,  can  be  written  as: 

u,  v,  w  =  [U  ( z),V  (z),W  (z)]cos(m4>  +  ^  )  Equation  3 

for  n> 0.  From  the  above  equation,  it  is  obvious  that  the  cylinder  supports  standing  waves  in 
the  circumferential  direction,  with  n  integer  waves. 

Excitation  of  the  structural  waveguide  can  be  accomplished  via  embedded  smart  material 
actuators.  The  field  of  embedded  smart  material  actuators  is  quite  extensive  .  Piezoelectric 
actuators  can  be  bonded  to  a  structure  and  can  produce  an  in-plane  strain  on  the  structure 
when  driven  electrically.  There  is  a  plethora  of  research  on  these  actuators  bonded  to  structures 
and  they  have  been  applied  to  such  diverse  fields  as  active  noise  control  and  control  of 
helicopter  blade  aerodynamics.  Excitation  of  a  cylinder  by  a  piezoelectric  actuator  is  derived  in 
Lalande  .  The  piezoelectric  actuator  is  viewed  as  a  distributed  in-plane  strain  forcing  function. 
To  arrive  at  a  particular  solution  for  a  cylinder  with  rigid  boundary  conditions,  this  forcing 
function  is  decomposed  into  the  modal  eigenfunctions  of  the  cylinder,  generating  modal 
participation  factors  for  the  actuator.  To  accomplish  a  propagating  traveling  wave  in  a  cylindrical 
structure  excited  by  embedded  smart  material  actuators,  phased  array  actuators  were 
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employed  .  The  main  concept  is  that  multiple  sensors/actuators  can  be  used  in  unison  to 
create  an  array  that  is  particularly  sensitive  in  one  spatial  domain,  in  this  instance,  the  modal 
domain.  To  excite  specific  structural  modes,  several  actuators  can  be  driven  in  a  spatial  phased 
array  .  These  actuators  will  be  driven  with  a  spatial  phase  array  concept  where  the  signal 
driving  one  actuator  will  be  delayed  by  a  time  (At)  that  will  reinforce  the  wave  propagating  from 
the  previous  actuator  .  By  repeating  this  concept  N  times,  the  phased  array  actuator  generates 
only  one  traveling  wave  corresponding  to  the  actuator  spacing  and  time  delay. 

Non-circular  structures  are  finite  and  will  therefore  exhibit  modal  behavior.  Recently,  there  have 
been  efforts  to  model  waveguides  with  arbitrary  cross-sections,  periodic  structures,  and  in- 
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homogeneities  .  Therefore,  one  significant  challenge  is  to  design  non-circular  structures 
which  posses  these  properties  e.g.  similar  shape  to  an  airfoil.  One  potential  solution  is  to  apply 
embedded  smart  materials  that  can  create  and  absorb  waves  in  the  finite  structure  that  will 
match  the  impedance  of  the  “infinite"  structure  at  the  boundaries.  This  can  also  be 


accomplished  passively  using  a  distributed  vibration  absorber3*,  which  can  effectively  absorb 
specific  wavenumbers  in  the  structure.  However,  an  investigation  into  non-cylindrical 
waveguides  was  not  pursued  because  of  the  time  and  financial  constraints  of  this  project. 


STRUCTURAL  WAVEGUIDE  AND  TEST  RIG  DESIGN 


The  test  rig  has  been  designed  in  order  to  represent  as  close  as  possible  a  simply  supported 
cylinder.  For  this  configuration,  an  analytical  solution  is  available  for  the  normal  modes. 

In  the  literature  it  is  possible  to  find  different  formulations  for  the  description  of  the  cylinder 
dynamics  that  are  valid  for  several  boundary  condition  and  cylinder  properties.  In  this  study  we 
have  focused  on  thin  walled  structures  where  the  ratio  between  the  thickness  of  the  cylinder 
and  its  diameter  is  small. 

If  we  consider  an  infinite  cylinder,  the  equation  of  motion  for  axial,  radial  and  torsional  direction, 
are  highly  coupled  unlike  the  case  of  plates  and  beams.  The  solution  in  terms  of  displacement 
in  the  three  principal  direction  is  assumed  to  be: 

u(x„  t )  =  Unm  cos(n)e**~‘**n=i  m  =1 
v(x,,t )  =  V„msin(n)e'**"'n=om=i 
w(x„t )  =  Wnm  cos(n)e**~n  =0  m  =1 

where  n  and  m  correspond  respectively  to  the  azimuthal  and  axial  modal  order.  Substituting  in 
the  equations  of  motion,  assuming  no  force  acting  on  the  system,  we  obtain  an  eigenvalue 
problem  whose  solution  is  usually  plotted  in  the  so-called  dispersion  curves  (wave  number  kns 
as  a  function  of  frequency). 

In  the  case  of  finite  length  cylinders,  the  procedure  previously  described  can  be  applied  to 
derive  modes  shapes,  natural  frequencies  and  forced  response  of  the  system.  In  this  case,  we 
consider  a  cylinder  of  length  L  that  is  simply  supported  at  the  two  ends  (radial  displacement, 
bending  moment,  torsional  displacement  and  shear  force  in  these  points  are  zero).  The  radial 
displacement  is: 


ifi 

w  (x,,t )  =  2  jW  cos(n)  sin (k  x  )e 
and  the  corresponding  eigenvalue  is: 
k  =  mO  /  L 


that  again  can  be  used  to  plot  the  dispersion  curves  for  this  configuration. 

This  analysis  is  fundamental  for  the  design  of  the  test  rig  in  order  to  evaluate  the  frequency 
range  of  interest  and  the  normal  modes  that  can  be  excited  using  the  piezoelectric  actuators. 
Moreover,  in  order  to  measure  any  changes  in  the  aerodynamic  drag,  it  is  necessary  to 
constrain  the  cylinder  in  such  a  way  that  it  is  possible  to  measure  the  displacement  in  the  axial 
direction.  This  can  be  realized  by  supporting  the  cylinder  on  tow  points  along  its  axis  and 
inserting  a  load  cell  in  between  the  constraint  points  and  the  supporting  structure. 

The  cylinder  chosen  for  this  work  has  the  following  properties: 


Low  thickness  to  diameter  ratio  (0.009). 

8”  diameter  to  guarantee  access  to  the  inside  and  enough  room  for  pressure  probes, 
accelerometers,  actuators  and  load  cells. 

Aluminum  alloy. 

These  characteristics  also  assure  that  the  circumferential  modes  that  we  are  interested  in  are  at 
low  frequency  as  confirmed  by  the  modal  hammer  analysis  and  the  plot  of  the  dispersion  curves 
for  the  cylinder  shown  in  Figure  1.  The  circumferential  mode  order  is  equivalent  to 
circumferential  wavenumber  and  is  calculated  by  performing  a  spatial  Fourier  transform  on  the 
measurements  around  the  cylinder.  The  red  sections  of  the  contour  plot  represent  regions  of 
high  level  and  correspond  to  circumferential  modes  i.e.  where  an  integer  number  of 
wavelengths  fit  into  the  circumference.  These  modes  can  then  be  targeted  as  potential  drive 
frequencies  for  the  traveling  wave  system. 
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Figure  1:  Measured  mode  order  (or  equivalently  circumferential  wavenumber)  vs  frequency  for 
_ the  cylinder  used. _ 


TESTING  AND  CHARACTERIZATION  OF  STRUCTURAL  WAVEGUIDE 
DESIGN 


This  document  describes  briefly  the  test  rig  used  to  demonstrate  the  possibility  to  drive 
circumferential  traveling  waves  in  a  cylinder  that  interacts  with  the  boundary  layer  in  order  to 
reduce  turbulent  drag.  After  a  description  of  the  test  rig,  the  theory  behind  generating 
circumferential  waves  in  cylindrical  structures  is  presented  and  the  system  identification 
procedure  is  illustrated. 

TEST  RIG  DESCRIPTION 

The  test  rig  is  an  aluminum  pipe  described  in  Figure  4.  The  cylinder  is  cantilevered  on  a  rigid 
support.  It  is  equipped  with  7  piezoceramic  actuators  equally  spaced  along  a  circle  at  a  location 
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that  is  half  of  its  length.  The  cylinder  is  internally  supported  by  means  of  4  load  cells  on  a  rigid 
beam  as  depicted  in  Figure  2.  The  load  cells  are  placed  in  such  a  way  to  measure  the  axial 
load.  In  order  to  uncouple  the  shear  drag  on  the  cylinder  from  the  pressure  acting  on  the  front 
surface,  the  front  blunt  body  has  been  mounted  directly  on  the  supporting  beam  by  means  of  a 
screw  and  the  gap  in  between  the  cylinder  and  the  blunt  body  has  been  filled  with  grease.  A 
detailed  drawing  of  the  front  part  of  the  test  rig  is  reported  in  Figure  3. 
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Figure  2.  Detail  of  load  cell 


Figure  3.  Detail  of  front  blunt  body 


Figure  4.  Test  rig:  lateral  section 

CIRCUMFERENTIAL  WAVES  IN  CYLINDERS 
Every  actuator  /  is  driven  with  a  sine  wave  and  produces  a  pair  of  waves  traveling  in  opposite 


direction  whose  expression  is: 


wi  (,t  )=  e 


+  +  jm*  +  yt  ♦  jm\ 


+ 


Equation  4 


^  ^  jm*  +  yt  +  jm\  i 

w  (,t  )=  e 


The  quantity  w  in  Equation  5  is  the  radial  displacement  and  is  a  complex  number;  in  order  to 
have  the  real  displacement  it’s  sufficient  to  take  its  real  part.  The  term  0  is  the  angular 
coordinate  with  respect  to  the  positive  x  axis,  w  is  the  angular  frequency  and  (pi  is  the  angular 
position  of  the  /- th  actuator.  The  parameter  m  is  the  order  of  the  circumferential  mode.  Then 
Equation  4  is  valid  only  for  normal  modes  since  the  wave  spatial  pattern  has  a  wavelength  that 
is  a  multiple  of  the  perimeter.  In  this  case  w  is  the  natural  angular  frequency  of  the  mode.  When 
the  system  is  excited  with  a  harmonic  function  at  co  the  two  waves  generate  a  standing  wave 
that  is  by  definition  a  normal  mode  of  the  system.  Each  actuator  is  able  to  excite  the  same 
mode  but  since  it  has  a  (pi  angle,  normal  modes  have  the  same  shape  (if  m  is  fixed)  but  they  are 
shifted  by  -mcpi.  In  the  most  general  case  of  N  actuators  driven  at  the  same  frequency  but  with 
arbitrary  phase  delay  yi the  expressions  for  the  two  waves  in  Equation  4  are: 


+  jm *  +  yt  ♦  ymi  *  i 


Equation  5 


*  ❖  jm ♦  +  yt  +  jmi  +  i 

Wj  (,t  )=  e 


And  when  are  all  driven  at  the  same  time,  the  resulting  displacement  takes  the  form: 


w  (.t  )=  w.\, t  )+w  *  (, t ) 


NN 


Equation  6  w  (*,t)=*e 
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Rearranging  Equation  6,  we  have: 


NN 


yt  +  jm 

Equation  7  w  (,t  )=  e  e  e 
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Figure  5.  Examples  of  circumferential  modes  (k  =  2, 3, 4, 5) 

TRAVELING  WAVES 

In  order  to  have  a  traveling  wave,  the  radial  displacement  in  every  point  on  the  circumference 
has  to  be: 


w  (,t  )=  ❖  cos (m^>  +  t ) 

this  can  be  seen  as  the  real  part  of  a  complex  signal  such  as: 

w  (,t  )=  Re(cos(m^  + 1 )  +  j  sin(mO  +  t )) 

that  in  the  complex  exponential  notation  is: 

w  (,t  )=  Re(e  *"*) 

This  is  equivalent  to  a  positive  traveling  wave  generated  by  an  actuator  placed  at  ip  =0.  Then,  in 
order  to  generate  a  traveling  wave,  we  have  to  suppress  the  one  of  the  two  waves  (depending 
on  the  desired  spinning  direction).  This  can  be  done  opportunely  phasing  the  signal  that  drives 
each  actuator;  this  means  that  the  first  sum  in  (4)  should  go  to  1  while  the  second  should  go  to 
zero.  Since  the  two  terms  are  equal,  except  for  their  signs,  we  choose  to  put  to  zero  every 
exponent.  In  this  way  the  two  sum  goes  to  N  when  we  set: 


Equation  8 

Then  equation  (3)  becomes: 


w  (,t  )=  Re(N  (cos (m)  +  j  sin(m)  +  cos(m)  ❖  j  sin(m))e  ) 


❖  i 


=  m\ 


Equation  9 


jm^ 


and  the  real  displacement  is: 


Equation  10  w  (,t  )=  Re(2N  cos(m)(cos(t )  +  j  sin(t ))  w  (,t  )=  2N  cos (m)  cos(t ) 


In  order  to  establish  the  real  properties  of  the  cylinder,  a  system  characterization  has  been 
performed  when  the  system  is  in  its  final  configuration  (cantilevered  on  rigid  support,  with  all  the 
sensors  and  wires  installed  and  the  front  and  rear  blunt  bodies).  The  scope  of  the  system 
identification  is  mainly  to  plot  the  dispersion  curves  and  verify  how  the  supports  are  influencing 
the  ideal  free-free  boundary  conditions.  Two  different  tests  have  been  carried  out: 

modal  roving  tap  test  (for  dispersion  curves) 

system  transfer  function  (point  acceleration  along  cylinder  axis  and  piezoelectric  actuator 


excitation) 


In  this  test  the  cylinder  transfer  functions  have  been  measured  in  a  point  on  the,  tapping  with  a 
modal  hammer  in  21  equally  spaced  points  along  the  circumference  of  the  cylinder  (this  is  the 
so-called  roving  technique  and  is  well  founded  under  the  assumption  that  the  principle  of 
reciprocity  is  valid).  The  force  and  acceleration  data  is  sampled  at  8192  Hz  and  the  transfer 
functions  are  averaged  over  5  tests  and  collected  (every  transfer  function  is  computed  using  a 
8192  coefficients  FFT  in  order  to  have  a  1  Hz  resolution  on  the  frequency  spectrum).  The 
dispersion  curves  are  derived  computing  the  spatial  FFT  on  the  collected  data. 
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Circumferential  modal  order 


Figure  6.  Dispersion  curve 


SYSTEM  TRANSFER  FUNCTION  ALONG  AXIS 


In  order  to  verify  the  system  response  along  the  axis  of  the  cylinder,  the  transfer  function 
between  the  acceleration  taken  in  24  positions  and  the  voltage  driving  one  of  the  piezoceramic 
actuators  have  been  taken.  The  results  are  presented  in  Figure  i.  Transfer  function  along  axis  where 
location  1  is  referred  to  the  cylinder  aft  while  24  is  the  front. 


Figure  7.  Transfer  function  along  axis 


It’s  clear  that  the  first  and  the  second  circumferential  modes,  respectively  at  129  and  355  Hz 
have  the  same  axial  wave  number  (equal  to  one)  while  other  modes,  for  example  at  500  Hz, 
have  an  higher  ka  that  result  in  nodal  lines  at  certain  axis  location.  Observing  the  system 
response  in  those  points  will  lead  to  a  poor  signal  while  position  12  is  a  good  observation  point 
for  all  modes. 


EXPERIMENTAL  SETUP 


Every  actuator  is  driven  by  a  PCB  790  series  amplifier.  Every  channel  can  drive  a  capacitive 
load  such  as  a  piezoceramic  actuator,  at  200  V  peak  to  peak  and  200  mA.  The  frequency 
bandwidth  depends  on  the  load  itself  and  with  the  installed  actuators  it’s  up  to  1.2  kHz.  The 
control  software  has  been  developed  in  C  on  a  Linux/RTAI  platform.  A  GTK+  user  interface 
controls  the  real  time  code  and  sets  up  the  communication  with  two  data  acquisition  boards 
(namely  a  National  Instrument  Nl  6143  for  analog  to  digital  conversion  and  a  Nl  6733  for  analog 
outputs).  The  code  allows  setting  the  sampling  frequency,  the  number  of  channels  to  drive  and 
the  parameters  of  the  driving  signal  (amplitude,  frequency  and  phase).  The  phase  is 
automatically  calculated  while  the  frequency  is  derived  from  a  configuration  file  that  has  been 
compiled  according  to  the  system  identification  procedure  (and  in  particular  from  the  dispersion 
curves  in  Figure  6  whose  summary  is  here  reported  in  table  1. 

Table  1.  Modes  frequency  and  order 


Mode  order 

Frequency  Hz 

2 

129 

3 

355 

4 

745 

PRELIMINARY  FLOW  VISUALIZATION  AND  PARTICLE  IMAGE 
VELOCIMETRY  RESULTS 


PARTICLE  IMAGE  VELOCIMETRY  BACKGROUND 

The  range  of  applications  and  current  advancements  of  DPIV  is  so  vast  and  diverse  that  a 
detailed  review  is  impossible  in  itself.  We  will  however  address  the  DPIV  contributions  that 
pertain  to  the  development  of  our  in-house  developed  TRDPIV  methodology  that  was  used 
during  this  effort.  Review  articles  by  Adrian  and  Grant  provide  comprehensive  summaries  of 
the  principles  of  the  method  and  its  various  applications.  The  work  of  Willert  and  Gharib  , 
Westerweef 41,  and  Huang  et  al.  established  the  digital  implementation  of  PIV.  Further 
developments  resulted  in  super-resolution,  iterative,  and  hybrid  DPIV  algorithms  that  increased 
accuracy  and  spatial  resolution  .  For  iterative  DPIV  methods  the  use  of  a  discrete-window 
offset  (DWO)  was  introduced  .  The  combination  of  iterative  DPIV  with  DWO  significantly 
reduces  the  spatial  averaging  effects,  improves  resolution  and  accuracy,  and  renders  the 
uncertainty  independent  of  the  total  measured  displacement.  The  accuracy  is  further  improved 
by  the  use  of  an  adaptive  Central  Difference  (CDI)  second  order  DWO  .  This  is  especially 
important  for  vortical  flows  and  systems  with  large  pulse  separation  such  as  pDPIV  and 
TRDPIV  systems. 


Despite  the  advances  in  improving  the  spatial  resolution  and  accuracy  of  the  method, 
conventional  DPIV  systems  provide  a  limited  temporal  resolution  (<30Hz),  which  is  insufficient 
for  resolving  turbulent  flow  fluctuations.  Previous  implementation  of  time-resolved  systems 
employed  digital  cameras  with  limited  resolution  or  drum  film  cameras  that  are  limited  in 
recording  times  .  Vlachos  presented  one  of  the  very  first  fully  digital  systems  that  extended 
the  sampling  frequency  to  1  kHz  with  256x256  pixel  spatial  resolution.  Recently  Upatnieks  et 
al.  presented  a  kHz  frame-rate  analog  PIV  system  capable  of  4000  frames/s  (fps)  for  0.5s 
maximum  recording.  Analog  recording  delivers  superior  spatial  resolution  (>1Kx1K).  However, 
due  to  the  high  frame  rates  the  extreme  rotation  speed  for  the  film  compromises  the  image 
registration.  This  significant  error  source  is  superimposed  on  the  common  digitization  errors. 
Vlachos  and  collaborators  recently  developed  a  TRDPIV  system  with  1-10  kHz  sampling  rate 
and  extended  recording  time,  capable  of  resolving  the  spatio-temporal  evolution  of  turbulent 
flow,  improving  upon  the  earlier  works  of  Scarano  &  Riethmuller  ,  Wereley  &  Meinhart  , 
Guezennec  &  Kiritsis  ,  and  Cowen  &  Monismith  .  Improvements  on  the  method  increased  the 
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accuracy  for  resolving  vortical,  near-wall  and  poly-dispersed  multi-phase  flows  .  We  introduced 
a  novel  hierarchical  scheme  that  is  computationally  efficient  and  has  showed  great  potential  for 
high  shear  vortical  flows  by  alleviating  the  large  spatial  averaging  effects  of  multigrid 
approaches.  A  new  hybrid  particle  tracking  method  based  on  an  adaptive  cross-correlation 
masking  was  also  introduced.  Detailed  error  analysis  by  our  group  characterized  the 
performance  of  the  system  and  compared  its  accuracy  with  commercial  and  other  established 
algorithms  .  Error  analysis  results,  using  Monte-Carlo  simulations  of  an  ideal  vortex,  comparing 
our  method  with  conventional  multigrid  DPIV  showed  that  the  overall  accuracy  especially  under 
high  shear  rates  was  improved  by  60%.  These  results  clearly  illustrated  the  advantages  of  our 
methodology,  particularly  for  resolving  vortical  flows. 

Although  during  the  past  few  years  commercial  TRDPIV  systems  have  become  available, 
successful  implementations  of  the  method  are  limited.  One  reason  is  the  prohibitive  cost  of 
these  TRDPIV  systems.  More  importantly,  TRDPIV  measurements  are  not  trivial.  Severe 
limitations,  including  the  very  low  pulse  energy  of  the  available  high  frequency  lasers  and 
accurate  timing  of  the  time  separation  between  image  pairs,  restrict  the  applicability  of  the 
method  to  a  few  types  of  -mainly  low  speed-flows  and  require  meticulous  experimentation.  It  is 
beyond  the  scope  of  this  text  to  provide  a  detailed  technical  description  of  these  issues.  The  PI 
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was  one  of  the  first  to  deliver  kHz  TRDPIV  measurements  over  a  variety  of  applications 
and  has  published  over  90  refereed  conference  and  journal  papers  related  method 
development  and  implementation. 

After  two  equipment  upgrades  we  are  currently  operating  a  stereoscopic  TRDPIV  system  with 
0-IOkHz  sampling  rate  and  Ips  minimum  inter-frame  time  delay.  In  addition,  3  IDT  X-Stream 
CMOS  cameras  with  1280x1024  pixel  resolution  at  1  kHz  sampling  rate  are  used  for  the  image 
recording.  CMOS  cameras  are  preferred  for  TRDPIV  systems  as  they  deliver  better 
performance  at  high  frame  rates  at  a  much  lower  cost  compared  to  CCD’s.  Moreover,  signal  to 
noise  ratio  of  current  CMOS  sensors  is  very  comparable  to  CCDs.  The  Stereo  DPIV 
implementation  used  in  our  system  is  based  on  the  work  by  Prasad  &  Adrian  .  In  conclusion, 
the  DPIV  system  employed  herein  is  one  of  very  few  available  worldwide  to  deliver  kHz  rate 
Time-Resolved  Stereo  DPIV,  it  involves  the  latest  available  instrumentation  and  it  extensively 
tested  and  benchmarked  over  a  wide  variety  of  flows. _ 


PRELIMINARY  TRDPIV  EXPERIMENTAL  VALIDATION 


Due  to  time  constrains  and  limitation  of  resources,  at  the  time  of  this  report  the  TRDPIV  study 
was  not  completed  therefore  all  results  shown  here  are  merely  representative  of  the 


measurement  capabilities.  Ongoing  experimental  results  and  further  analysis  of  these  efforts 
will  be  disseminated  at  a  later  time  in  the  form  of  conference  proceedings  and  archival 
publications. 


EXPERIMENTAL  METHOD 


A  TRDPIV  experimental  study  was  performed  in  order  to  quantitatively  visualize  the  effect  on 
the  flow  of  driving  circumferential  travelling  waves  around  the  wall  of  a  cylinder.  The  effect  of 
travelling  waves  under  quiescent  flow  conditions  was  investigated.  The  experimental  methods 
and  procedures  are  explained  and  some  preliminary  findings  are  demonstrated. 

An  in-house  customized  TRDPIV  system  was  employed  in  order  to  quantitatively  visualize  the 
flow  very  close  to  the  wall  of  the  travelling  wave  cylinder.  A  series  of  instantaneous  velocity 
fields  were  generated  using  this  technique.  Figure  8  (left)  provides  a  view  of  the  experimental 
setup.  It  includes  a  New  Wave  Nd-YAG  laser,  which  produces  approximately  1.5  centimeter 
wide  and  2  millimeter  deep  laser  sheet  with  the  use  of  properly  placed  optics.  The  cylinder  is 
located  inside  a  plenum  in  order  to  keep  the  particle  seeding  dense  enough  for  PIV.  Flow  was 
seeded  with  pulverized  90%  solution  of  Di-2-Ethylhexel  Sebecat  using  the  Lavision  seeder 
shown  in  Figure  8  (right). 


Figure  8.  (a)  View  of  the  TRDPIV  experimental  setup  (b)  cylinder  structure  and  seeder 


Implementing  the  method  described  earlier,  a  structural  travelling  wave  is  applied  around  the 
circumference  of  an  aluminum  cylinder  in  a  clockwise  direction  by  the  use  of  7  equally  spaced 
piezoceramic  actuators.  The  cylinder  was  forced  with  a  travelling  wave  at  its  third 
circumferential  mode  which  is  demonstrated  in  Figure  9. 
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Figure  9.  The  third  circumferential  mode  shown  above  was  used 

An  IDT  MotionPro  X5  digital  CCD  camera  was  used  to  acquire  images.  A  special  optical 
component  called  a  GRINscope  was  used  with  the  camera.  This  is  an  in  house  developed  and 
customized  borescopic  imaging  system  based  on  a  cylindrical  lens  with  radially  varying  index  of 
refraction.  This  enables  it  to  transmit  the  light  with  high  efficiency  and  minimal  distortion.  This 
component  made  it  possible  to  achieve  a  magnification  of  approximately  30  pm  per  pixel 
without  needing  to  place  the  camera  close  to  the  laser  sheet.  A  closer  view  of  the  GRINscope 
and  the  region  of  interest  is  shown  in  Figure  10.  The  resulting  magnification  was  sufficient  to 
resolve  the  micron  size  structural  deflections  on  the  surface  of  the  cylinder. 


RESULTS 


Figure  11  through  Figure  13  show  the  velocity  fields  at  three  successive  time  instants. 
At  the  first  instant,  shown  in  Figure  1 1 ,  the  cylinder  wall  is  forced  by  the  travelling  wave 
to  move  to  the  location  designated  by  the  dashed  curve.  As  expected  there  is  a  strong 
flow  away  from  the  wall.  In  the  next  time  instant  shown  in  Figure  12,  the  cylinder  wall 
retreats  to  its  original  location  designated  by  the  solid  curve,  and  there  is  a  flow  in  the 
clockwise  direction  around  the  circumference  of  the  cylinder  wall.  At  the  time  instant 
shown  in  Figure  13,  the  flow  decelerates  and  reaches  also  quiescent  conditions. 


Figure  1 1 .  First  time  instant,  fluid  is  being  pushed  away  the  wall  by  the  wall’s 

movement. 

Images  were  acquired  at  200Hz  sampling  rate.  Acquired  images  were  processed  by  using  the 
FlowlQ  software  developed  by  Aeroprobe  cooperation.  FFT  based  correlation  technique  was 
employed  using  16x16  pixel  interrogation  windows.  Iterative  and  discrete  window  offset  (DWO) 
method  were  used  along  with  adaptive  central  difference  DWO  algorithms  in  order  achieve 
higher  accuracy  for  a  near-wall  flow  condition. 

RESULTS 

Figure  11  through  Figure  13  show  the  velocity  fields  at  three  successive  time  instants.  At  the 
first  instant,  shown  in  Figure  1 1 ,  the  cylinder  wall  is  forced  by  the  travelling  wave  to  move  to  the 
location  designated  by  the  dashed  curve.  As  expected  there  is  a  strong  flow  away  from  the  wall. 
In  the  next  time  instant  shown  in  Figure  12,  the  cylinder  wall  retreats  to  its  original  location 
designated  by  the  solid  curve,  and  there  is  a  flow  in  the  clockwise  direction  around  the 
circumference  of  the  cylinder  wall.  At  the  time  instant  shown  in  Figure  13,  the  flow  decelerates 
and  reaches  also  quiescent  conditions. 


Figure  13.  Third  time  instant,  flow  is  steady. 


Figure  12.  Second  time  instant, there  is  a  flow  in  the  clockwise  direction. 
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CONCLUSIONS  AND  FUTURE  WORK 


During  the  nine  month  performance  of  the  work  presented  herein  we  were  successful  in 
designing  developing  and  fabricating  structural  waveguides  within  a  cylindrical  structure. 
Structural  characterization  of  the  prototype  demonstrates  that  the  performance  was  within  the 
design  expectation.  Subsequently  aerodynamic  testing  was  initiated  in  order  to  validate  that  the 
proposed  concept  can  impart  flow  disturbances  that  reduce  turbulent  viscous  drag.  At  the  time 
of  this  report  only  quantitative  flow  visualization  under  quiescent  flow  conditions  had  been 
performed  and  preliminary  analysis  of  the  results  is  presented.  The  PIV  results  showed  a  three 
stage  cycle  occurring,  which  creates  a  motion  that  sweeps  the  flow  very  close  to  the  wall  in  a 
clockwise  direction.  This  movement  of  fluid  might  interrupt  the  near-wall  turbulence  production 


and  reduce  the  viscous  drag.  Future  work  on  this  project  should  include  drag  measurements  in 
an  air  tunnel  in  order  to  see  the  effect  of  travelling  wave  on  the  drag  measurements. 
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